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ABSTRACT 

We developed a quantum field theoretical description for the surface states of 
three-dimensional topological insulators. Within the relativistic quantum field the¬ 
ory formulation, we investigated the dynamics of low-lying surface states in an 
applied transverse magnetic field. We argued that, by taking into account quantum 
fluctuations, in three-dimensional topological insulators there is dynamical gener¬ 
ation of a gap by a rearrangement of the Dirac sea. By comparing with available 
experimental data we found that our theoretical results allowed a consistent and 
coherent description of the Landau level spectrum of the surface low-lying excita¬ 
tions. Finally, we showed that the recently detected zero-Hall plateau at the charge 
neutral point could be accounted for by chiral edge states residing at the magnetic 
domain boundaries between the top and bottom surfaces of the three-dimensional 
topological insulator. 
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1 Introduction 


Topological insulators realize new quantum states which have recently attracted consider¬ 
ably interest in condensed matter physics (for recent reviews, see Refs. [H |2l [3l IH [5l |6l [TJ 
[8],[9l[T0] and references therein). Topological insulators can be realized in both two dimen¬ 
sions and three dimensions. In particular, topological insulators in three dimensions m 
are non magnetic materials having an energy band gap in the bulk and possessing metal¬ 
lic surface states due to the nontrivial topology of the bulk electronic wavefunctions. In 
this paper we focus on time-reversal invariant systems, where the nontrivial topology is 
protected by time-reversal symmetry. In fact, in three-dimensional topological insulators 
the electron spin is locked to the momentum due to time-reversal symmetry leading to 
the notion of helical surface states. In these systems it turns out that the surface states 
display Dirac dispersions. Thus, the physics of low-lying excitations are described by 
relativistic Dirac fermions. As a result, the surface electronic band structure is similar to 
that of graphene (see, for instance. Refs. [niiiaiiiiiig), except that there is just a single 
Dirac fermion instead of four as in graphene (due to valley and spin degeneracies). Ap¬ 
parently it seems that the description of the two-dimensional surface of three-dimensional 
topological insulators in terms of a single Dirac fermion violates the Nielsen-Ninomiya no- 
go theorem [16] which states that for a time-reversal invariant system Dirac points must 
always come in pairs. In fact, the partner Dirac fermion resides on the opposite surface 
of the solid, restoring the counting scheme for Dirac fermions. Finally, the surface states 
of a topological insulator cannot be localized even for strong disorder as long as the bulk 
band gap remains intact. As a result, the surface states are topological protect against 
disorder and non-magnetic impurities. Remarkably, the existence of these surface states 
has been recently conhrmed by angle-resolved photoemission spectroscopy and scanning 
tunneling microscopy [D S H EHZ). 

The low-energy excitations in three-dimensional topological insulators are described by an 
effective Hamiltonian made of a two-component Pauli spinors which satisfy the massless 
two-dimensional Dirac equation with the speed of light replaced by the Fermi velocity 
vp- If the topological insulator is immersed in a transverse magnetic held, the relativistic 
massless dispersion of the electronic wave functions results in non-equidistant Landau 
levels (cgs units): 

= Ed ± ^ 2 n h — eB , n = 0,1,2,... (1.1) 

where eB > 0, e being the elementary charge, and Ed is the energy of the charge-neutral 
Dirac point. Note that, here and in the following, we shall not distinguish between the 
magnetic hux density B and the magnetic held H since in a nonmagnetic material B = H. 
We see that the Landau quantization of massless Dirac fermions is characterized by the 
occurrence of Landau levels, called zero modes, pinned at the Dirac point with energy 
Eq = Ed and the appearance of states with energies scaling with y/n on both the pos¬ 
itive and negative energy sides of the Dirac point. In fact, the presence of anomalous 
Landau levels at the Dirac point leads to the half-integer quantum Hall ehect [l2]. Re¬ 
markably, quite recently such peculiar Landau quantization of the surface states in three- 
dimensional topological insulators has been conhrmed by means of the scanning tunneling 
spectroscopy. 
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The plan of the paper is as follows. In Sect. |2] describe the dynamics of the low-lying 
surface states by means of an effective relativistic quantum held theory. In Sect. [3] we 
discuss the quantum dynamics of the surface states in presence of a transverse magnetic 
held. Sect. IHis devoted to the discussion of the dynamical generation of a mass gap by 
a rearrangement of the Dirac sea. In Sect. [5] we compare our theoretical results with 
available experimental data. In Sect. | 6 ]we briehy discuss the quantum Hall ehect, and 
suggest that chiral edge states could explain the zero-Hall plateau recently detected in 
three-dimensional topological insulators. Finally, our conclusions are relegated in Sect. [71 
For reader’s convenience some technical details are collected in Appendix [^ 

2 Quantum Field Theory of Topological Surface States 

In this section we are interested in the surface states of three-dimensional topological 
insulators. More formally, we shall describe the dynamics of these surface states by means 
of an ehective relativistic quantum held theory, namely quantum electrodynamics in two 
spatial dimensions. It turns out that the held theory description of such states should be 
adequate for our purposes. We, further, assume that the chemical potential is pinned at 
the charge-neutral Dirac point. Without loss in generality we may set /i = Ejj = 0. Thus, 
the low-lying excitations are described by the following ehective Hamiltonian [Dial Em 

[61 E]: 

H = ( 2 . 1 ) 

where a = 7^7 and vp is the Fermi velocity. Here 7 ^ = ( 7 °, 7 ^, 7 ^) are the Dirac gamma 
matrices which satisfy the Clihord algebra : 

{7^7^} = 2^/^% (2.2) 

= diag{l, —1, —1) being the Minkowski metric tensor. In two spatial dimensions a 
spinor representation is provided by two-component Dirac spinors. Then, the fundamental 
representation of the Clihord algebra is given by 2 x 2 matrices which can be constructed 
from the Pauli matrices as follows: 


7° = 7^ 

= icTi, 7 ^ = ia2 . 

(2.3) 

The Heisenberg equations of motion are: 



. ^ d'^{x,t) 

m 

= , 

(2.4) 

. ^ (f,f) 

dt 

= ¥{x,t),H . 

(2.5) 


From the equations of motion together with the canonical equal time fermion anticom¬ 
mutation relations: 


|T(a;, t), f)| = 6{x — x') , |T(a;, t), T(a;', f)| = t), f)| = 0 , 

( 2 . 6 ) 
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one can determines the quantum dynamics of the system. In fact, we write the held 
operators in terms of the creation and annihilation operators (see, for instance, the classic 
textbook Ref. HZl): 


T (T, t) = 

J d^p^^e * « bp -1- « 'ijj^ \x) dt| , 

(2.7) 

= J 

f bp + 4 } > 

(2.8) 


where are the positive and negative energy solutions of the well known Dirac 

equation (see Appendix and the creation and annihilation operators satisfy the anti¬ 
commutation relations: 

V'} = > {4’ 4 -'} = ^iP-P’) (2-9) 

with all the other anticommutators vanishing. The bt and bp operators create and destroy 

particle above the Fermi sea, while dt and dp operators create and destroy hole inside 
the Fermi sea. We rewrite the Hamiltonian in terms of the the creation and annihilation 
operators. To this end we insert Eqs. fl2.7p and fl2.8p into Eq. fl2.ip . Using Eqs. flA.4D and 
flA.Sp we obtain: 

H = j d?pep ^bVbp - dp 4} = j d^P {^> h + 44} + ^0 , (2.10) 

where: 

Eo = - j epd^p-p) . (2.11) 

To properly interpret the d-function in Eq. fl2.1ip . we note that: 

SO that: 

^ / (2^ ^ (2^’ 

where V is the two-dimensional volume (area) of the system. Accordingly we have: 

The Hamiltonian Eq. fl2.10p is now normal ordered with respect to hlled Fermi-Dirac sea. 
In fact. Eg given by Eq. fl2.14p is the energy of the ground state (the vacuum) dehned by: 

6p-|0 > = dp-|0 > = 0 • (2-15) 

It is now evident from Eq. fl2.15p that the ground state corresponds to the hlled Fermi- 
Dirac sea. If we try to evaluate the vacuum energy Eq, we face with the problem of 
ultraviolet divergencies. In fact, using ep= vp |p| we obtain: 
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As is well known, the problem of divergent qnantnm corrections is inherent to relativistic 
quantum held theories. The inhnities encountered in relativistic quantum theories can 
be cured with a procedure called renormalization d- Actually, in condensed matter 
physics one never had to deal with ultraviolet divergences. In fact, the description of the 
low-lying excitations by means of the effective Hamiltonian Eq. fl2.10p is valid as long as 
ep < W, where W is the band-width. Therefore there is a natural ultraviolet cut-off 
which makes hnite the vacuum energy Eq. (12.161) : 


Eo ~ 


- 1 / 


1 

Gvr fPvp 


(2.17) 


This (negative) vacuum energy contributes to the binding energy of the whole system 
and, therefore, does not inhuence the dynamics of the low-lying excitations. In this case 
the renormalization corresponds to simply subtract Eq from the Hamiltonian operator H 
leading to the renormalized Hamiltonian: 


Hr = H - Eq 


(fp Ep 


bij 
P P 


+ 


dp dp 


(2.18) 


This renormalization corresponds to the so-called normal ordering prescription. For the 
purposes of the present work we do not have to do any further renormalization. 

For later convenience, we rewrite the ground-state energy as: 


Eq = - V 


d^p 

(2nfiy 



(2.19) 


To evaluate the integral we use the following computational trick which has been has been 
already employed by us [T8] : 


We get: 


Eq = 




ds d 


TIS ds 


e 


—as 


d^p r ± _yj: _ ^ r ds d n\ 

(27rh)^ Jq ■\/7rs ds dvr fp Vp Jq ^/7Fs ds \s J 


( 2 . 20 ) 


( 2 . 21 ) 


Equation fl2.2ip shows that the ultraviolet divergency manifest itself as a singularity for 
s —?• 0. Introducing an effective band-width W such that s > l/hE^, we easily obtain: 


Eq ~ 


- V 


1 

3 7lyEr fpvp 


( 2 . 22 ) 


Comparing this last equation with Eq. fl2.17p we see that the cut-off W is hnitely related 
to the cut-off W,W = 


3 Topological Fermions in Magnetic Fields 


A distinguish feature of two-dimensional Dirac fermions is the peculiar Landau quantiza¬ 
tion of the energy states in applied magnetic helds. In presence of an external magnetic 
held the effective Hamiltonian operator Eq. fl2.ip becomes: 


H = 


d^x T^(a;, t) Ivr a ■ 


—ih'V -I- - A{x) 
c 


T(x,f) , 


(3.1) 
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where the electromagnetic vector potential is such that V x A{x) = B{x). To write the 
Hamiltonian operator in terms of the creation and annihilation operators, it is convenient 
to expand the fermion held operator f) and T(T, t) in terms of the wave functions 
which are the solutions of the Dirac equation in presence of the external magnetic held. 
In other words, we shall adopt the so-called Furry picture |T9]. Thus, we expand the 
fermion held operator Tl(T, t) and T(T, t) in terms of the wave function basis ijj^p and 
given in Appendix 


r’+co 


, n=l 


T(f,t) = / Kp + dip + V’o,p(^) A 


p ( 1 

(3.2) 


f + CxD 


. n=l 


4't(f,t)=/ V’iftU *>l,p + « ■"* 4,p +i/'J,p(S)4 


(3.3) 

The creation and annihilation operators operators satisfy the standard anticommutation 
relations: 


. Kp] = |5(P - p') p, , ipp] = |5(P - p') Sp.n' , [il, , Cpj = S(p - p') , 

(3.4) 

all the other anticommutators vanishing. Inserting Eqs. fl3.2p and fl3.3p into the Hamilto¬ 
nian operator Eqs. fl3.1|) we hnd: 


^ p-\-oo ^ ^ ^ ^ ^ r-\-oo ^ ^ ^ ^ 

H = ^ ^ / dpSn I dn,p ^ ^ / dp£n | + *^n,p J TEo(i?) , 

n=l d-oo 2-00 ^ 

(3.5) 


with: 


°° P + OO TD °° 

Eo{B) = -^ dpend{p-p)dn,n = - V 
n=l 2-00 


(3.6) 


The physical interpretation of Eq{B) is straightforward. In fact, from Eq. fl3.5p we infer 
that the ground state of the Hamiltonian is given by: 


bn,p \0;B> = dn,p \0;B > = 0 , 


(3.7) 


SO that: 

H \0-, B > = Eo{B) \0-, B > (3.8) 

namely, Eq{B) is the vacuum energy in presence of the external magnetic held B. 

Since for large n ~ ^/n (see Appendix |^, we see that the summation over n in Eq. fl3.6p 
is divergent. In other words, Eq{B) is ahected by ultraviolet divergences. However, if we 
consider the renormalized Hamiltonian Eq. fl2.18p we hnd: 

'+00 

dpSn 

where: 

AEo(H) = Eo{B) - Eo . (3.10) 



CO 

Hr = H - Eo = Y 
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Now we show that AEq{B) is not affected by any ultraviolet divergences. In fact, in 
Appendix 1X1 we show that: 


AEo{B) 


Eo{B) — Eq 


eB 

2Trhc 


ds d 

1 

1 

/o y/Esds 

1 - 

2hvl^s 

r c 


(3.11) 


This last equation shows that, indeed, AEo{B) is finite. After some elementary manipu¬ 
lations we rewrite Eq. flS.lip as: 


cB I cB 

AEoiB) =^VvF^2n— X /, (3.12) 

where: 

j _ f°° dx d 
Jq y/nxdx 

We end, thus, with the remarkable result that in presence of an applied magnetic field 
the ground state due to quantum fluctuations acquires a finite positive energy given 
by Eq. fl3.12p . To understand the physical meaning of AEo{B), we observe that the 
degeneracy of the Landau levels, namely the number of states in a given level, is: 


1 — e“ 


X 


dx 


(1 - e-^)^ 


1 

+ —17 


X^ 


~ 0.208 . (3.13) 


9l 


eB 

2tiHc 


V . 


Therefore, we may define the vacuum energy per particle as: 


(3.14) 


e^(B) = ~ 0.208 VF \kn— . (3.15) 

9l V c 

Now, it is quite easy to convince ourself that Sp{B) modifies the single particle spectrum 
Eq. fll.ll) as: 


= Ef) -I- Sp{B) ± y2fiVF — n , n = 0,1,2, ... (3.16) 

Usually one assumes that the external magnetic field does not modify the Dirac point en¬ 
ergy. However, the subtle quantum effects we have evaluated show that external magnetic 
fields give rise to a shift of the Dirac point. Remarkably, this effect manifest itself also in 
graphene [211122]. We believe that the effect has been already observed in the experimen¬ 
tal investigation of the quantum Hall effect. In fact, it turns out that by switching on an 
external magnetic field one must tune the gate voltage Vg to drive the system to the Dirac 
neutral point. Moreover, it seems that the gate voltage needed to reach the neutral point 
depends on the magnetic field as Vg ^/H, in qualitative agreement with our Eq. (I3.15p . 

4 Dynamical Generation of Mass Gap 

In Sect.[2]we said that in three-dimensional topological insulators the low-lying excitations 
are described by the effective Hamiltonian Eq. (12.ip . Indeed, this Hamiltonian corresponds 
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to massless two-dimensional Weyl fermions. In two dimensions we can define the parity 
and time-reversal transformation [2S1 UM [2S] according to: 


where x 


{xi,X 2 ) and x' 


VA\x,t)V-^ 
VA\xA)V-^ 
VA^{x,t)V^^ 
V V-^ 

(-Xi,X 2 ), 

ta\xa)t-^ 
T A {x,t)T-^ 
T^{x,t) r-^ 


= A%x',t) 

= —A^{x',t) 

= A^{x'A) 

= ai-^iyX'A) , 

= A^{x, —t) 

= —A{x, —t) 

= a2'^{x,-t) 


(4.1) 


(4.2) 


It is, now, easy to show that the Hamiltonian Eq. fl2.ip is invariant under V and T 
transformations. However, if one is interested in massive fermions, then the relevant 
Hamiltonian turns out to be: 


H 


(fx ^^(x, t) i—i hvF a -V + mv%a^\i>{xA) ■ 


(4.3) 


It is easy to recognize that the mass term is odd under both V and T transformations. 
It is important to point out that this parity-violating mass is, in fact, the only possibility 
for two-dimensional Weyl fermions. 

Long time ago we showed [IB] that in (2-1-1 )-dimensional quantum electrodynamics with 
two-component Dirac fermions in an external magnetic helds it is energetically favored to 
develop a fermion mass gap. In the present Section we shall show that this happens also 
to the surface states of three-dimensional topological insulators. 

In presence of applied magnetic helds the relevant Hamiltonian is given by Eq. fl3.ip . 
Since we are interested in massive fermions, then the effective Hamiltonian becomes: 


H= (fx'^^{x,t)\vF<y- 


—ih'V + - A{x) -I- mvpa^ [>\l/(a;, f). (4.4) 


The results of Ref. [18] would imply that it is energetically favored to develop a gap such 
that: 

m Vp = Aq sign{B) , (4.5) 

with Aq > 0. Actually, the magnetic held is a vector. However, we are considering 
magnetic helds which are perpendicular to the surfaces of the specimen. Thus, as con¬ 
cern the dynamics of the surface quasiparticles the magnetic held can be considered a 
(pseudo)scalar. Our convention is that R > 0 if the magnetic held is directed along the 
normal to the surface, while R < 0 otherwise. Note that under parity and time-reversal 
transformations both R and m change sign, so that the gap Aq in Eq. (14. 5 p is V and T 
invariant . 

To be dehnite, we hx the direction of the magnetic held such that R > 0. Accordingly we 
are left with the following Hamiltonian: 


H= (fx'^\x,t)^VF<y- —ihV + - A{x) + Aq as | T(T, f) . (4.6) 
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To write the Hamiltonian in the Furry representation, we need the solutions of the Dirac 
equation: 


Vf o. ■ 


-itW + -A{x) + A 0 CT 3 = s'tpix) . 


Proceeding as in Appendix we get: 


( 4 . 7 ) 


I Aq 
~ Vp 

-hdi + i—X2 


hdi — i 


,-e_B 

C 


ihdo 


X2 - ihd2 

Vp 


iIj{x) = — iIj{x) 
Vf 


Using Eq. flA.9p we get: 


[ip - i^X2 - ind2\ <j>2(X2) = 4 ^ 1 ^ 1 ( 12 ) 
[-ip + i^X2 - iW2] lt>l{x2) = ^^'t>2(x2) ■ 


( 4 . 8 ) 


( 4 . 9 ) 


From Eq. fl4.9p we easily obtain: 


-fa? + 


1 6=^2 


2 - 1^2 W 2“^(^2 - 01 ( 3 ^ 2 ) - I 

h{x2) = [-ip + ipX2 - ihd2] (j)l{x2) . 


£2-A2 




heB 

c 


(f>l{x2) 


(4.10) 


The solutions of Eq. fl4.10p are given in Eq. flA.lTp with eigenvalues: 
£ = ib 


eB 


2h-—{n + l)u^ + Aq , u = 0, 1, 


(4.11) 


Obviously, we have also the zero modes with wavefunctions given by Eq. (lA.lSp and 
eigenvalue: 

£0 = - Ao . (4.12) 

Therefore, we are lead to the following spectrum: 




(+) 



2 h—+ Aq 
2h—+ Aq 


n 


n 


1 , 2 ,... 

0 , 1 , ... 


(4.13) 


where the eigenfunctions are given by Eq. flA.21F For convenience, in Fig. [T]we compare 
the spectrum of the Landau levels with A = 0 to the case of non-zero gap. 

Thus, we can write: 


T (T, t) 




OQ 



n =0 



( 4 . 14 ) 



. n=l 




gn t 
h 


.0(+)t| 

T n,p ' 


x) b^p + * 


,• gn i 


n=0 



( 4 . 15 ) 


where the creation and annihilation operators operators satisfy the standard anticommu¬ 
tation relations, and 


gn 


2h— VpU Al 


(4.16) 
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Figure 1: Schematic spectrum of Landau levels of surface states in topological insulators 
in applied magnetic field (left). Landau levels with dynamical generation of the gap A 
(right). 


Whereupon, we obtain the Hamiltonian operator: 


H 



71=1 


h 

^n,p 


71=0 


d 

^n,p ^n,p 


+ E,{B) 


with: 


Eo{B) = - V 


eB 


CO 


27rhc ^ 

71=0 


The renormalized Hamiltonian operator is: 


( 4 . 17 ) 


( 4 . 18 ) 


Hr = 


/ +CO [ ^ 

I n=l 


b 

^n,p ^n,p 


+ 


71=0 


d 

^n,p ^ri,p 


+ AEo{B) 


(4.19) 


where now: 

AEo{B) = - H ^ ^ ^2n^vln + Ag - E^ . (4.20) 

n =0 

To evaluate the (divergent) sum in Eq. fl4.20p we employ the same trick as before. After 
some algebra it is easy to show that: 


CO 


E 

n =0 



+ Ao 


ds d 

1 - 

(MO 

<1 

1 

1 _ 

Jo V^s ds 

^+2nvl^s _ ^ 


(4.21) 
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Figure 2: Plot of the function g{\), Eq. fl4.27p . versus A. 


Since 


we readily obtain: 


e-Ag. 


s^O 
~ + 


AEo{B) = V 


eB 

2tiHc 


-An + 


ds d 
IFs ds 


2hvp^s 


e-Ag. 


^ 2nvl^s 


Therefore the vacuum energy per particle is: 
ep{B,Ao) = — Ao 


7°° ds d 

e-Ags 

1 

Jq ^/^TS ds 

^+2h.vj,^s _ ^ 

2nvi^s 


(4.22) 


(4.23) 


(4.24) 


Introducing the dimensionless variable: 


A = 



Eq. (I4.24j) can be rewritten as: 


ep{B, Aq) 



where: 


gW = -A + 


1'°° dx d 

H 

1 

1 _ 

1 

Jq \fT^ dx 

e* — 1 

X 


(4.25) 


(4.26) 


(4.27) 
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It is straightforward to see that the function g{\) is hnite. Indeed, in Fig. |2]we display 
g{X) versus A. Eqs. fl4.26p and fl4.27p imply that the vacuum energy per particle depends 
not only on the magnetic held but also on the gap Aq, which until now has been a free 
positive parameter. It is worthwhile to stress that for Aq = 0 Eq. fl4.26p agrees with 
Eq. fl3.15p since 5f(0) ~ 0.208. Looking at Fig. [2] we note that the vacuum energy per 
particle decreases for small values of the gap Aq, reaches a minimum, and then increases 
for large enough values of the gap. In fact, we hnd the minimum at A = A ~ 0.56 with 
g{X) ~ —0.093. To understand this peculiar behavior of the vacuum energy, we observe 
that the opening of the gap Aq affects the spectrum of the Landau levels in two way. On 
the one hand the gap decreases the energy of the zero modes (see Fig. [1]). On the other 
hand the energy of higher Landau levels tend to increase with the gap. It is the balancing 
of these two effects that results in the vacuum energy minimum. Thus, we see that it is 
energetically favorable to induce a dynamical gap by a rearrangement of the Dirac sea 
triggered by quantum fluctuations. This remarkable phenomenon is quite analogous to 
the Peierls |26] quantum phase transition in quasi one-dimensional conductors. 

In summary, we found that the low-lying surface states in three-dimensional topological 
insulators immersed in an external magnetic held give rise to gapped Landau levels: 


■(+) 

n 

— Ej:){B) + -1 

\j2h—Vpn + Aq , 

n = 1, 2, 

.(-) 

n 

= Ejd^B) — -1 

\j2h^Vpn + Aq , 

n = 0, 1, 


(4.28) 


with: 


Aq 


= A 



~ 0.56 



(4.29) 


and 

I eB 

Ed(B) = Eo ep{B , Aq) ~ Ep — 0.093 y 2Hvp — . (4.30) 

The results obtained in this Section rely on our effective Hamiltonian which should de¬ 
scribe the low-lying excitations of the surface states in topological insulators. One must 
still face the task of comparing the results obtained within a highly-idealized model with 
experimental investigations of three-dimensional topological insulators. 


5 Landau Levels: Comparison with Experimental Data 

In the preceding Section we argued that it is energetically favored by quantum huctu- 
ations to develop a mass gap for the low-lying surface excitations in three-dimensional 
topological insulator in an external magnetic held. The aim of this Section is to contrast 
our theoretical expectations to experimental data. 

Quite recently, several experimental studies conhrmed that the peculiar relativistic na¬ 
ture of surface states in topological insulators in magnetic helds leads to Landau levels 
which varies linearly with y/riB. The authors of Ref. [27l [2H] reported direct observation 
of Landau quantization of the topological surface states in Bi 2 Se 3 in magnetic helds by 
using scanning tunneling microscopy and spectroscopy. Relativistic Landau levels are also 
observed in the tunneling spectra in a magnetic held [29] by using scanning tunneling spec¬ 
troscopy on the surface of the topological insulator. Moreover, microwave spectroscopy 
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has been applied to study cyclotron resonance due to intraband transitions between Lan¬ 
dau levels in Bi2Se3 |3U] . 

The hallmark of Dirac fermions for low-lying surface excitations in three-dimensional topo¬ 
logical insulator is the presence of a held independent Landau level at the Dirac point 
(zero mode) and the already alluded scaling of higher Landau levels with \/nB. In fact, 
Fig. 4 in Ref. 127 ] displays the energies of the Landau levels versus \/nB for magnetic 
helds ranging from 8 T to 11 t[§. Remarkably, for higher values of n one observes Landau 
level energies which vary linearly with \/ nB. Moreover, from Fig. 2 in Ref. one infers 
that the lowest Landau level, assumed to be the zero mode n = 0, is pinned at the Dirac 
point: 

Ed ~ —200 mev , (5.1) 

consistently with the value inferred from the minimum of the differential tunneling con¬ 
ductance in absence of magnetic held. However, for low n, namely for the Landau levels 
with energies near the Dirac point, it seems that the Landau levels deviate from the linear 
behavior. We believe that this ehect can be accounted for if we admit the presence of a 
dynamical mass gap according to our previous discussion in Sect. 01 To check quantita¬ 
tively this, we need the energies of the Landau levels for diherent magnetic held strength. 
Fortunately, Fig. 2 in Ref. [28] reports the Landau levels for diherent magnetic helds 
ranging from 6 T up to 11 T. We used the experimental data extracted from Fig. 2 of 
Ref. [2S] to perform a quantitative comparison with the theoretical results discussed in 
the previous Section. To this end, we note that in Refs. [23 l2H] the energy of the Lan¬ 
dau levels are determined by htting the diherential conductance with multiple gaussians. 
The typical peak width of Landau levels turns out to be of about 8 mev. Therefore, to 
be conservative, for our analysis we assumed for the Landau level energies a statistical 
uncertainty of about 2-3 mev. 

Preliminarily, to check our procedure, we htted the data to the expected spectrum of 
Landau levels for massless Dirac fermions, Eq. jnii. Since the average Fermi velocity in 
Bi2Se3 lies in the range vp = 3.0 — 7.0 x 10^ cm/s [31], we write vp = x 5.0 10^ cm/s. 
Accordingly, for the spectrum lying above the Dirac point, Eq. (HU can be rewritten as: 

= Ep + 18.34 mev \/ n B{T) , ?7, = 0,1,2,.., (5.2) 

where B(T) stands for magnetic held strength in Tesla. 

To implement the hts we used the program Minuit [32] which is conceived as a tool to hnd 
the minimum value of a multi-parameter function and analyze the shape of the function 
around the minimum. The principal application is to compute the best-ht parameter 
values and uncertainties by minimizing the total chi-square 

We htted the experimental data to Eq. fl5.2p for each values of the magnetic held leaving 
Ep and as free htting parameters. Moreover, to get a sensible ht to the data we ex¬ 
cluded the Landau levels with n = 0,1. In Tabled] we summarize the results of the htting 
procedure. Note that in the last column of Table dl we indicate the reduced chi-square, 
namely the total chi-square divided by the number of degree of freedom. As rule of 
thumb, a sensible ht results in xl ~ 1- 

In Fig. |3]we display the experimental data together with the htting lines for two rep¬ 
resentative values of the magnetic held. As it turns out, the lowest Landau level has an 

^Even though we are using cgs units, it is widespread practice to measure the strength of the magnetic 
field in Tesla. We recall that 1 T = 10^ G. 
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Figure 3: Landau level energies e^\B) versus \/nB for B = 7 T (Top) and B = 11 T 
(Bottom). Data have been borrowed from Fig.2 of Ref. [28]. The assumed statistical 
errors are smaller than the symbol size. Dashed lines are the hts of e^\B) to Eq. fl5.2p . 
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Table 1: Summary of the values of fitting parameters in Eq. (15.211 reported, respectively, 
in the second and third columns. In the fourth column we give the reduced chi-square. 


B(T) 

Ed ( 

mev) 

V5 



Xr 

6 

-231 

± 

5 

1.27 

± 

0.03 

1.20 

7 

-237 

± 

5 

1.29 

± 

0.03 

1.00 

8 

-249 

± 

5 

1.32 

± 

0.03 

1.05 

9 

-246 

± 

5 

1.30 

± 

0.03 

0.77 

10 

-251 

± 

5 

1.32 

± 

0.03 

0.95 

11 

-256 

± 

6 

1.36 

± 

0.03 

0.99 


energy quite close to the Dirac point Ed ~ - 200 mev independently on the magnetic held 
strengths. However, the lowest Landau levels n = 0,1 seems to deviate appreciably from 
the linear behavior ~ \/nB which hts nicely the n >2 Landau levels. Even more, if we 
extrapolate e^n'^ according to Eq. fl5.2p we obtain the energies of the Dirac point, showed 
in Table dl which, not only lie very far from the value at zero magnetic held, Eq. fl5.ip . 
but seem to decrease with the magnetic held. On the other hand, the htted values of the 
Fermi velocity vp does not show a sizable dependence on the magnetic held strengths. In 
fact, in Fig. |l]we report the htted values of the Fermi velocity versus the magnetic held 
strength. We hnd for the average Fermi velocity: 

vf = 6.55 ± 0.12 10^ cm/s , xl = 1-07 . (5.3) 

It is reassuring to see that nj? in Eq. fl5.3p lies in allowed range for the average Fermi 
velocity in Bi 2 Se 3 [37] . 

The puzzling anomalies encountered within the standard interpretation of the measured 
Landau level energies point to an alternative interpretation of the experimental data. We 
shall, now, show that the measurements can be accounted for in a consistent manner 
within the theoretical scenario discussed in Sect. 01 In this case, the spectrum of the 
Landau levels is given by Eqs. fl4.28p . fl4.29p . and fl4.3Up . To this end, we need the Landau 
levels which are above the Dirac point. According to Eq. fl4.28p we have: 


with: 


and 


For convenience we rewrite Eq. fl5.4l) as: 


— vj,n + Ag , n = 1, 2, ... , 

( 5 . 4 ) 

'1— , S ~ 0.56 , 
c 

( 5 . 5 ) 

/ eB 

y 2 fiVp — , a ~ 0.093. 

( 5 . 6 ) 

^/B{T) [n + 52) , n = 1, 2, ... . 

( 5 . 7 ) 


We have htted the available data to Eq. (15.7p leaving Ed{B), v^, and 6 as free parameters. 
As a result, the htting procedure gives the best ht values of the parameter for each given 
value of the magnetic held. For convenience, the results of our hts are summarize in 
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Figure 4: Fermi velocity in 10^ cm/s versus the magnetic field B(T). The dashed line is 
the average velocity, Eq. fl5.3p . obtained by fitting with a constant the values in Tabled] 


Table 12] while in Fig. [5] we display the experimental data together with the fitting curves 
for three representative values of the magnetic field. 

A few comments are in order. First, we note that the Dirac energy Ed{B) do depend 
on the magnetic field as in the previous fits. However, presently we expect a shift of the 
Dirac point according to Eq. fl5.6p . In fact, as we will discuss later the shift of the Dirac 
point due to the magnetic field seems to be consistent with theoretical expectations. From 
Table |2] we see that the Fermi velocity is independent on the magnetic field strengths. In 
fact, we find for the average Fermi velocity: 

vf = 7.2 ± 0.1 10^ cm/s , = 1-67 . (5.8) 

This value of the Fermi velocity is slightly higher, but in reasonable agreement with our 
previous determination Eq. fl5.3p . 


Table 2: Summary of the values of the fitting parameters in Eq. fl5.7p reported, respec¬ 
tively, in the second, third and fourth columns. In the last column we give the reduced 
chi-square. 


B(T) 

Ed{B) (mev) 

V5 



h 



xi 

6 

-262.0 

± 

4.2 

1.39 

± 

0.03 

0.50 

± 

0.36 

1.33 

7 

-271.1 

± 

4.4 

1.42 

± 

0.02 

0.50 

± 

0.36 

1.21 

8 

-283.5 

± 

4.4 

1.44 

± 

0.03 

0.50 

± 

0.42 

0.92 

9 

-288.9 

± 

4.5 

1.45 

± 

0.02 

0.50 

± 

0.35 

0.67 

10 

-298.1 

± 

4.6 

1.48 

± 

0.02 

0.50 

± 

0.35 

0.76 

11 

-299.3 

± 

4.9 

1.49 

± 

0.08 

0.55 

± 

0.13 

0.42 
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[nB(T)]''' 


Figure 5: Landau level energies Sn versus \/nB for B = 6 T (Top), B = 7 T (Middle), 
and B = 11 T (Bottom). Data have been extracted form Fig.2 of Ref. [28]. The assumed 
statistical errors are smaller than the symbol size. Dashed lines are the fits of e\^\B) to 
Eq. dEH). 
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Figure 6: Dirac point energies E^, in meV versus the magnetic field B in Tesla. Dashed 
line is the fit of Ed{B) to Eq. fl5.6l) . 


As concern the spectrum of the Landau levels, we stress that the levels lying above the 
Dirac point have n > 1, according to Eq. (15.7p . In fact, the zero mode n = 0 has been 
pushed below Ed^B) due to the dynamical generation of the mass gap. Therefore, Fig. [5] 
shows that our theoretical curves are able to track all the experimental data. Note that 
the lowest Landau level n = 1 has an energy which seems to be almost independent on 
the magnetic field. This is due to the compensation of two different effects. Indeed, the 
increase of the energy due to the magnetic field for n 7 ^ 0 is contrasted by the negative shift 
of the Dirac point. For n = 1 these two effects almost perfectly compensate each other. 
From Fig. |5] and Table [2] we may conclude that the data are in reasonable agreement with 
Eq. fl5.7p for all the available values of the magnetic field. Note that in Fig. [5l to magnify 
the non linearity in \/nB due to the gap, the fitting curves have been draw starting from 
nB = 0. Concerning the parameter <5, looking at Table [2] we see that this parameter does 
not depend on the magnetic field. We find for the average value: 

5 = 0.51 ± 0.18 , xl = 0.01 , (5.9) 

which is in remarkable agreement with our theoretical expectations. 

Finally, we must check if the dependence of the Dirac energy on the magnetic field is 
in agreement with Eq. fl5.6p . In Fig. | 6 ] we report the fitted values of the Dirac energy 
in Table [2] as a function of the magnetic field. We have fitted these values according to 
Eq. (15.op assuming for the Fermi velocity the average value Eq. (15.8p . and leaving Eo and 
a as free parameters. The fit returns for these two parameters the following values: 

Ed = - 190 ± 14 mev , a = 1.22 ± 0.04 , xl = 1-53 . (5.10) 
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Indeed, we find that the shift of the Dirac point with the magnetic field is in satisfying 
agreement with Eq. fl5.6p (see dashed line in Fig. E]). Moreover, the parameter Eq, which 
is the energy of the Dirac point in absence of the magnetic field, is in good agreement 
with the experimental value Eq. fIS.ip . However, we must stress that the parameter a is 
well above our theoretical estimate, Eq. f|5.6p . We recall that this parameter quantify the 
vacuum energy per particle. In our opinion, since the calculations in previous Section have 
been performed within a highly idealized model, such a discrepancy does not spoil the 
overall coherent agreement of our theoretical picture with the experimental observations. 


6 Quantum Hall Effect and Chiral Edge States 


The Hall effect indicates the voltage drop across a conductor transverse to the direction 
of the applied electrical current in the presence of a perpendicular magnetic field. In the 
quantum Hall effect, with increasing magnetic field, the Hall resistance evolves from a 
straight line into step-like behaviors with well-defined plateaux. At the plateaux the Hall 
conductance is quantized 

( 6 . 1 ) 


^xy — 


2tiH 


where u is an integer or a certain fraction. At the same time, the longitudinal resistance 
drops to zero, suggesting dissipationless transport of charged quasiparticles. 

A single two-dimensional Dirac fermion under a magnetic field is known to show the 
quantized Hall effect with the Hall conductance: 


^xy 


& 1 

M + 2> ■ 


(6.2) 


n being an integer. The factor 1/2 in Eq. fl6.2p is characteristic of the relativistic dispersion 
relation of Dirac fermions compared with the usual massive electrons, and it is related to 
the zero modes pinned at the neutral point. Alternatively, the half-integer quantization 
can be also accounted for as a Berry phase [3Sl El]- The Hall conductance of three 
dimensional topological insulators is expected to be given by the sum of the contributions 
from the top and bottom surfaces (see Fig. [T] left panel). Therefore one finds: 


(^xy = {n + m + 1) , (6.3) 

27rh 

with n and m are integers referring to the upper and lower surface respectively. When 
the Landau levels of the top and bottom surface coincide exactly, the two contributions 
are equivalent and one has n = m, namely only the odd integer quantum Hall effect is 
expected. If, however, the degeneracy of the Landau levels is removed, then the Hall 
conductance becomes quantized as integers. Indeed, recently the integer quantum Hall 
effect has been observed in three dimensional topological insulators [35l 136] . 

If we admit the generation of a dynamical mass gap, then, according to Eq. (14. 5 h we obtain 
the gapped Landau level spectrum for the two surfaces as displayed in Fig. [T] right panel. 
In fact, we adopted the convention that magnetic fields perpendicular to the surfaces are 
positive in the direction of the outward surface normal, so that the transverse magnetic 

^ In this Section we shall indicate the spatial coordinates as {x,y) instead of {xi,X 2 )- 
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Figure 7: Schematic spectrum of surface state Landau levels in applied magnetic field 
(left). Landau levels with dynamical generation of a gap Aq (right). The Fermi level is 
at Dirac neutral point = 0. 


held is positive for the top surface and negative for the bottom surface. It is evident that 
also in this case we obtain the odd integer quantum Hall effect for identical surfaces, or 
the integer quantum Hall effect if the the degeneracy of the Landau levels in the top and 
bottom surfaces is removed. 

Interestingly enough, in high magnetic helds an additional quantum Hall plateau at the 
charge neutral point with u = 0 has been observed [36]. This peculiar zero-Hall plateau 
could be explained in a natural way if there are gapless edge excitations pinned at the 
Dirac neutral point. Indeed, it has been already proposed [37] that chiral edge states 
residing at the magnetic domain boundaries are responsible for the zero-Hall plateau 
observed in the quantum anomalous Hall effect in magnetic topological insulators. More 
recently, chiral edge states residing at the magnetic domain boundaries on the surface 
of a three-dimensional topological insulator have been observed by means of a scanning 
superconducting quantum interference device [38] . 

In this Section we discuss the presence of current-carrying chiral edge states residing 
at the magnetic domain boundaries between the top and bottom surfaces of a three- 
dimensional insulator. Moreover, these peculiar edge excitations are pinned at the Dirac 
neutral point, whereas the energy spectrum of the surface states is gapped. It turns out 
that these states are localized on the magnetic domain wall and are similar to the Jackiw- 
Rebbi soliton [33] . 

To face this task, in Fig. [8] we show a schematic picture of the magnetic held on the 
upper and lower surfaces of a three-dimensional topological insulator. At the edges of the 
surfaces the magnetic held decreases to zero within a region of linear size given by the 
magnetic length: 
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Figure 8: Schematic picture of the magnetic kink relevant for the edge states. Dots and 
crosses indicate the direction of magnetic held leaving and entering the surface, respec¬ 
tively. The width of the strip (dashed lines) where the magnetic held change sign is set 


by the magnetic length qb 



The physical meaning of Eq. (16.41) is as follows. The low-lying quantum states localized 
on the upper surface will feel a magnetic held which is positive and vanishes near the 
surface boundaries within a distance of order ~ a^. Obviously, the same happens on the 
lower surface where the magnetic held is negative. If we glue together the two surfaces 
along an edge, as shown in Fig. El we see that the surface excitations will be subject to 
an ehective magnetic held which can be written as: 


B{x) 


Bo tgh 



(6.5) 


where Bq is the strength of the (constant) magnetic held inside the two surfaces. Accord¬ 
ingly, the potential vector is given by: 

Mx,y) = {-B{x)y , 0) . (6.6) 


Let us consider, now, the Dirac equation in presence of the magnetic held given by 
Eq. fl6.5p . After taking into account Eq. fl6.6p . we readily obtain: 


-Hd^ + i^B{x)y 


ihd„ 


f^x — i^B{x)y — ihdy 
0 


i^{x,y) = —i){x,y), (6.7) 

Vp 
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Note that in Eq. fl6.7p we do not consider the mass gap Aq. In fact, we are interested here 
in the zero-energy solutions that are localized at the surface boundary x = 0 where the 
magnetic held vanishes. Let Q{x,y) be the solutions of Eq. fl6.7p with £ = 0. Writing: 


Q{x,y) 


0i{x,y) \ 

02{x,y) ) ’ 


( 6 . 8 ) 


we get: 


+fid^ - i^B{x)y - iHdy\ 92{x,y) = 0 
-Hd^ + i^B{x)y -ihdy ei{x,y) = 0 . 


Let us consider, hrstly, the region x > 0. For x S> evidently B{x) ~ Bq. In this case 
we have already seen that the solution of Eq. (16.91) is given by Eq. (I3.14p . This suggests 
to seek the solutions of Eq. (16. 9 p in the form: 


= ( 9(x,!,) ) • ^ > 0 • (6-10) 

Inserting Eq. (I6.10p into Eq. (16.9p we obtain the equation for 6{x,y) which we rewrite as: 

( 6 . 11 ) 

c; 

To solve this last equation, we write: 

e{x,y) = o{x,y) , (6.12) 


+fid^ - i^[B{x) - Bo\y - i^Boy - iHdy^ 9{x,y) = 0 . 


to get: 


b + CxD 


[B{x') — Bo\dx' — ihdy > 9{x,y) = 0 


+tld^ - i-Boy - 

c c 


To proceed further, we put: 

9{x,y) = 9p{x,y) = /(x) <I'o,p(|/) , 

where: 


^o,p( 2 /) = e , C = 

After some manipulations, we hnally get: 

fix) 


eBp 

he 


y 


cp 

eBn 


fix) 


. sW = -'"[1 + . 


giving: 


/(x) = A/ eV nc Jo V Ha . 

To a good approximation, for x < a^, we can rewrite Eq. (I6.17p as: 


f(x) ~ A/" e 


(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 

(6.18) 
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Putting it all together, we end with: 


Q+{x,y) = 


0 

Gp{x,y) 


, X > 0 , 


where: 

and 


It is easy to check that for a; < 0 the solutions of Eq. (I6.9jl are given by : 

0 -( 


e^(x,y) = _ 


e.(x,y) = ( ^ < 0 


(6.19) 

( 6 . 20 ) 
( 6 . 21 ) 


( 6 . 22 ) 


It is worthwhile to stress that the two zero-energy modes 0±(a;, y) are localized in the re¬ 
gion |x| < ^and are related by parity and time-reversal transformations. Therefore, 
if the chemical potential is tuned to the Dirac neutral point, then the zero-energy edge 
modes will contribute to the conductivity. In particular, it is evident that the contribu¬ 
tions of these modes cancel exactly in the Hall conductance a^y = 0, while they add to 
the longitudinal conductance a^x = CTyy > 0. So that, we see that these edge modes give 
rise to the Hall plateau at the charge neutral point with u = 0, but with a hnite longitu¬ 
dinal resistance at variance of the usual Hall plateaux where the longitudinal resistance 
vanishes. 


7 Summary and Conclusions 

In this work, we have discussed the dynamics of low-lying surface excitation in three- 
dimensional topological insulators. We developed a quantum held theoretical description 
for the surface states of three-dimensional topological insulators, which allowed us to 
investigate the quantum dynamics of low-lying surface states in presence of an applied 
transverse magnetic held. 

We evaluated the ehects of quantum huctuations on the ground state. In particular, we 
argued that, in presence of a constant transverse magnetic held, the quantum huctuations 
induce a shift of the energy of the Dirac neutral point which, in turns, depends on the 
strength of the magnetic held. Moreover, we argued that low-lying surface excitations 
in three-dimensional topological insulators develop a mass gap varying with \/eB by a 
rearrangement of the Dirac sea induced by quantum huctuations. Interestingly enough, 
very recently the physical consequences of dynamical mass generation in several classes 
of topological Dirac metals have been discussed in Ref. [ID] . 

To compare our theoretical scenario with observations, we reanalyzed the available ex¬ 
perimental data for the Landau level spectrum of the surface states in three-dimensional 
topological insulators. Remarkably, we argued that our theoretical results allowed a con¬ 
sistent and coherent description of the Landau level spectrum of the surface low-lying 
excitations. 

Finally, we showed that recently detected zero-Hall plateau at the charge neutral point 
could be accounted for by chiral edge states residing at the magnetic domain boundaries 
between the top and bottom surfaces of three-dimensional topological insulators. 
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A Appendix 

We collect here some intermediate steps needed in the derivation of results presented in 
Sects. |2]and[3l Firstly, let us consider the positive and negative energy solutions of the 
following Dirac equation: 


— ihvp a ■'V 

It is quite easy to solve Eq. flA.ip : 

1 


(x) = 


2tiH 


Up , V’J- \x) = ^ e Vf , ep = vf , 


where: 


Up — 


x/2 


-i 


, Up 


= ±( 


, tan. 9p = — . 


Pa 


(A.l) 


(A.2) 


(A.3) 


y/2 V +* 

Note that the positive and negative energy solutions of the Dirac equation are normalized 
according to: 

f (fx [ijj^^\x)]A^^\^) = ^{P~P')- (A.4) 


p' ' p 

Moreover, the Pauli spinors up and vp satisfy the following relations: 

1 0 


upui + vpvi = 1 = 


0 1 


, uiv_p = 0 . 


(A.5) 


With the aid of Eqs. flA.4p and flA.Sp one can verify that the anticommutation relations 
Eq. fl2.9p are indeed equivalent to the canonical equal time anticommutation relations 
Eq. (|2.6p . 

Next, we are interested in the Dirac equation in presence of an external magnetic held: 


Vp a ■ 


-ihV + -A(x) ipix) = e ipix) . 
c J 


(A.6) 


Since we are interested in uniform magnetic helds perpendicular to the crystal surfaces, 
in the Landau gauge we can write: 


Ak{xi,X2) = -X2 B 6k,i k = l,2. 


(A.7) 


As it is well known, in this case the Dirac equation Eq. flA.hp can be exactly solved (see, 
for instance. Ref. ESI). Indeed, Eq. (IA.6p leads to: 


0 hdi — — ifid2 

—Hdi + i^X2 — ihd2 0 


'0(F) = —'0(F) , (A.8) 

Vf 


with di = = 1,2. Writing: 


we obtain: 


)' 

(A.9) 

[ip - 02 ( 3 ^ 2 ) = ^01 ( 3 ^ 2 ) 

[-ip + - ih0a] 0i(a:a) = ^ 02 ( 2 ^ 2 ) • 

(A.IO) 
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Inserting the second eqnation into the hrst we rewrite Eq. flA.lOj) as: 

+ \^[ x , - mx 2 ) = I 01(0:2) 


02(0:2) = V [~^P + ^^^2 - ifid2] 01(0:2) . 


(A.ll) 


The hrst eqnation in Eq. (lA.lip is the familiar harmonic oscillator eqnation. Tlins, we 
can write: 


where: 


N„ = 


eB 


TrhcJ y'2^n\ 


snch that: 


^+00 


Moreover, the energy eingenvalnes are: 


eB 


One, also, easily hnd that: 

02 (^ 2 ) = ± ^n+l,p{^2) ■ 

Therefore, we are lead to the following eigenfnnctions: 


Nn HniO e--2^" , 

(A. 12) 

eB / cp \ 

Iw 0 ^ ~ 08 / ’ 

(A.13) 

• n. Note that onr 

normalization is 

= 1 . 

(A.14) 

n = 0, 1, ... 

(A.15) 


0S(a^na:2) = 


VM V2 V ±i^n+l,±pix2) 


(A.16) 


(A.17) 


with energy eigenvalnes given by Eq. flA.151) . respectively. It is easy to see that there also 
zero mode solntions: 

gipxi / g 

00,p (0:1, 0:2) = 


with energy eigenvalne: 


y/M V '^o,p(o:2) 

£0 = 0. 


If we adopt the convention that: 

*hn,p(o:2) =0 if n < 0 , 

then we can write: 

' VM V2 f ± »>!>„,±5(3:2) 

g+ipxi / g 

00,p (0:1,0:2) = 

with energy eigenvalnes: 


(A.18) 

(A.19) 

(A.20) 


n = 1, 2, 


V *^o,p(o:2) 


n = 0 . 


(A.21) 


£0 = 0 , = ±en = ±vf \l2H—n , n = 1, 2, 


eB 


(A.22) 
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The wave functions in Eq. flA.21D are normalized as: 


^ + 00 


= S{p-p') 6n,n' ■ 


(A.23) 


Note that the Landau levels are inhnitely degenerate with density of states: 

r-+oo , 


dp'lp^rSh^) = 


27ific 


(A.24) 


To derive Eq. fld.lip we note that by using Eq. (I2.20p we may write: 

eB r ds d 




— 


27iHc 




n=0 n=0 

We may, now, perform the summation over n: 


ns ds 


ns 


(A.25) 


n=0 


— 2hv%, — ns 


to obtain: 




1 — e 

ds d 


y/ns ds \i — 


(A.26) 


(A.27) 


As discussed in Sect. [21 the ultraviolet divergences in Eq{B) are recovered from the sin¬ 
gularity for s —0. To isolate the divergent part in Eq. (IA.27p . we note that: 


1 — e 

This suggests to rewrite Eq. flA.27p as: 

eB 


s^>0 , 1 

~ + 


2nvi^i 


(A.28) 


Eo{B) = 


271 he 
eB 


V 


ds d 

1 

1 

Jq y/ns ds 

\ _ g-2to|.^s 

2hvl^s 


27ihc 


1/ 


ds d 


ns ds 2hvp—s 


(A.29) 


The divergent term in Eq{B) is due to the second integral in the right-hand part of 
Eq. flA.29p . Comparing with Eq. fl2.2ip we at once recognize that this divergent term 
coincides with Eq. Therefore, we obtain: 


AEq{B) = Eq{B) -Eq = 


eB 

2nhc 




ds d 

1 

1 

Jq y/ns ds 

X _ 

2hvy-fs 


, (A.30) 


which agrees with Eq. fl3.1ip . 
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